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Abstract 


We consider the Schrddinger-Poisson-Newton equations for crystals with a cubic lattice and one ion per cell. 
We linearize this dynamics at the ground state and introduce a novel class of the ion charge densities which 
provide the stability of the linearized dynamics. This is the first result on linear stability for crystals. 

Our key result is the energy positivity for the Bloch generators of the linearized dynamics under a Wiener- 
type condition on the ion charge density. We also assume an additional condition which cancels the negative 
contribution caused by electrostatic instability. 

The proof of the energy positivity relies on a novel factorization of the corresponding Hamilton functional. 
We show that the energy positivity can fail if the additional condition breaks down while the Wiener condition 
holds. 

The Bloch generators are nonselfadjoint (and even nonsymmetric) Hamilton operators. We diagonalize 
these generators using our theory of spectral resolution of the Hamilton operators with positive definite energy 
021 [la. Using this spectral resolution, we establish the stability of the linearized crystal dynamics. 
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1 Introduction 


First mathematical results on stability of matter were obtained by Dyson and Lenard in lH |9l where the energy 
bound from below has been established. The thermodynamic limit for the Coulomb systems was studied first 
by Lebowitz and Lieb HU [191, see the survey and further development in ||20l. These results were extended 
by Catto, L. Lions, Le Bris and others to Thomas-Fermie and Hartree-Fock models ||4l|5l|6l. All these results 
concern either the convergence of the ground state of finite particle systems in the thermodynamic limit or the 
existence of the ground state for infinite particle systems. The dynamical stability of infinite particle ground 
states was never considered previously. 

We establish for the first time the dynamical stability of crystal ground state in linear approximation for the 
simplest Schrddinger-Poisson model. The ground state for this model was constructed in ifT^ . 

We consider crystals with the cubic lattice F = and with one ion per cell. The electron cloud is described 
by one-particle Schrddinger equation. The ions are described as classical particles that corresponds to the Born 
and Oppenheimer approximation. The ions interact with the electron cloud via the scalar potential, which is a 
solution to the corresponding Poisson equation. 

This model does not respect the Pauli exclusion principle for electrons. However, it provides a convenient 
framework to introduce suitable functional tools, which might be useful for physically more realistic models 
(Thomas-Fermie, Hartree-Fock, and second quantized models). In particular, we find a novel Wiener-fype 
sfabilify criferion (11.231) . (11.241) . 

This invesfigafion is mofivafed by fhe lack of a suifable mafhemafical model for a rigorous analysis of funda- 
menfal quanfum phenomena in fhe solid slate physics: heal conduclivily, eleclric conduclivify, Ihermoeleclronic 
emission, phofoeleclric effecl, Compton effecl, elc., see ill- 

We denote by a{x) fhe charge densify of one ion: 

[ o(x)dx = eZ>0, (1.1) 

where e > 0 is fhe elemenfary charge. Lef V^(x,f) be fhe wave function of fhe elecfron field, and <P(x) be fhe 
eleclrosfalic potential generated by fhe ions and elections. We assume h = c = m = I, where c is fhe speed of 


lighl and m is fhe election mass. Then fhe coupled equations read 




i\j/{x,t) = —-Ar/r(x,t) — e<F(x,f)t/r(x,f), 

X G 


(1.2) 

—A4>(v,t) = p{x,t) ■='^(y{x — n — q{n,t)) - 

n 


X G 

(1.3) 

Mq{n,t) = —{V^{x,t),(y{x — n — q{n,t))), 

n G Z^. 


(1.4) 


Here fhe brackels sland for fhe Hermilian scalar producl in fhe Hilbert space and for ifs differenl ex- 

lensions, and fhe series (11.31) converges in a suifable sense. All derivatives here and below are underslood in 
fhe sense of disfribufions. These equations can be written as fhe Hamilton system wifh a formal Hamilfon 
funclional ^ 

= i /J|Vt//'(v)|^ + p(v)Gp(v)]r/v + ^^^, (1.5) 

2 Jr3 “ 2M 

where G := —A^' and q := {q{n) : n G 7?), p := (p(n) : n G 7?), and p(x) is defined similarly to (11.31) . Namely, 
fhe sysfem (1 1 .2I) - (1 1 .41) can be formally written as 

/t/r(x,f) 1 qipd^ ^p(n)d^ 1 ^q(n)d^ 1 (1-6) 
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where dy := \{d^^ +*^ 2 ) with z\ = Rez and zi = Imz. A ground state of a crystal is a F-periodic stationary 
solution 

, <I>°(;c) , q^{n)=q^ for n£l? (1.7) 

with a real ft)° (and q^ € can be chosen arbitrary). A ground state was constructed in ifldll . Substituting (11.71) 
into (fOb-lTU]). we obtain the system 


xeT^:=M.^/r, ( 1 . 8 ) 

-A<I>°(;c) = p^{x):=a^{x)-e\\/'{x)\^, xeT\ (1.9) 

0 = — (V<I>°(;r),a(.r —« —^°)), n e (1-10) 

where we denote 

0^{x) :='^o{x — n — q^). ( 1 - 11 ) 

n 

In present paper we prove the stability for the formal linearization of the nonlinear system (ll.2b - (ll.4l) at the 
ground state (11.71) . Namely, substituting 

\ff{x,t) = [\ff°{x)+^'{x,t)]e^'‘^°‘, q{n,t) =q^ + Q{n,t) (1.12) 


into the nonlinear equations (11.21) . (11.41) with = Gp{xf), v/e.formally obtain the linearized equations 

(see Appendix A) 

[idt + ft)°]*F(x,t) = —j/^^{xf) — e<I>°(x)*F(x,t) — ei/r'’(x)Gpi {xf) 


Q{n,t) =P{n,t)/M 


X e 
n£l? 


(1.13) 


P{n,t) = — {VGpi{t), a{x — n — q^)) + ,'Va{x — n — q^)Q{n,t)) 


Here Pi{xf) is the linearized charge density 

Pi{x,t) = — ^ Va(x — « — q^)Q{n,t) — 2eRe [i/r'^(x)'R(x,f)], (1-14) 

n 

The system (11.131) is linear over M but it is not complex linear. This is due to the last term in (11.141) . which 
appears from the linearization of the term \ y\^ = ¥¥ <11-31) . However, we need the complex linearity for the 

application of the spectral theory. This why we will consider below the complexification of the system (11.131) 
writing it in the variables *Ti(x,t) := Re*T(x,t),*T2(-^T) •= Ini*T(x,t). We will consider the case when the 
ground state ¥^{x) can be taken to be a real function. In this case 


Re [r/r‘’(x)*T(x,t)] = r/r‘’(x)'Ti (x,t). 


Further we denote 

T(0 = ('Ri(-,0,'I^2(-T),<2(-T),^(-T))- 

Then (11.131 ) can be written as 


Y{t)=AY{t), 


/ 0 77 ° 

-77° - 2eV‘’G'/^” 0 

0 0 

\ -2S* 0 


0 0 \ 

-S 0 

0 

-r 0 / 


(1.15) 

(1.16) 


(1.17) 
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where := — — e^^{x) — (O^, the operators S and T correspond to matrices (14.41 ) and (14.51 ) respectively, 

and t//* denotes the operators of multiplication by the real function The Hamilton representation (11.61) 

implies that 


A = JB, 




/ + 0 

0 2H^ 

2S* 0 

\ 0 0 


2S 0 

0 0 

T 0 

0 M-^ 


\ 

/ 


(1.18) 


where J is the skew-symmetric matrix (15.21) . Our basic result is the stability for the linearized system (11.171) : 
for any finite energy initial state there exists a unique global solution, and it is bounded in the energy norm. 

We show that the generator A is densely defined in fhe Hilberf space ^ := L^(M^) 0L^(M^) 
and commufes wifh franslafions by vectors from F. Hence, fhe equafion (11.171) can be reduced by fhe Fourier- 
Bloch-Gelfand-Zak Iransform fo equafions wifh fhe corresponding Bloch generafors A(0) = JB{d), which 
depend on fhe parameter 6 from fhe Brillouin zone H* := [0,2?!]^. The Bloch energy operator B{9) is given by 


B{d) 


/ 2^0(0) 04eV°<5(0X 0 25(0) 0 \ 

0 2 ^°( 0 ) 0 0 

25*(0) 0 f(0) 0 

\ 0 0 0 / 


0 en*\r*, 


(1.19) 


where F* := 2kI?, and H^{d) := —^(V 0/0)^ — e<4>°(x) — ft)°. Further, G(0) is fhe inverse to fhe operator 
(/V — 0)^ : —)■ L^{T^). Finally, 5(0) and f{d) = ^ 2 ( 0 ) 0 fi (0) are defined respecfively by (17.221) and 

(Ejs, dm. 

However, the operator A is not selfadjoint and even not symmetric, which is a typical situation for the 
linearization of f7(l)-invariant nonlinear equations ifTSl Appendix B]. Respectively, the Bloch generators A(0) 
are not selfadjoint in the Hilbert space 


:=L2(^3)0L2(^3)0C^0C^ :=mVF. 


( 1 . 20 ) 


The main crux here is that we cannot apply the von Neumann spectral theorem to the nonself adjoint generators 
A and A(0). We solve this problem by applying our spectral theory of the Hamilton operators with positive 
energy imiisi, which is an infinite-dimensional version of some Gohberg and Krein ideas from the theory of 
parametric resonance ifT^ Chap. VI]. This is why we need the positivity of the energy operator B{d): 

^(0,F):=(F,B(0)F)r3>x(0)||F||^(^3), a.e. 0Gn*\r, (1.21) 

where >r(0) >0, the brackets stand for the scalar product in and we denote 

r(r^) :=H\T^)(SH\T^)eC^®C\ (1.22) 


This positivity allows us to construct the spectral resolution of A(0) which implies the stability for the linearized 
dynamics (11.171) . 

The key result of the present paper is the proof of the positivity (11.211 ) for the ions’s charge densities a 
satisfying the following conditions on the corresponding Fourier transform d{^). The first one is the Wiener- 
type condition 


Wiener Condition: 




L 




- ^=27Cm—9 


>0 


a.e. 0Gn*\F*. (1.23) 


3 























This condition is an analog of Fermi Golden Rule for crystals. The second condition reads 

a(27rm)=0, mGZ^\0. (1-24) 

The proof of the positivity (11.211) relies on a novel factorization of the Hamilton functional. This positivity 
necessarily breaks down at 6 G T*. Examples 18.11 and 18.21 demonstrate that the positivity can break down at 
some other points and submanifolds of IT*. 

Our main novelties are the following: 

I. The factorization of energy (16.41 ). (16.61 ) and (18.81 ). (18.101 ). 

II. The energy bound from below (16.11) for general densities o{x). 

III. The energy positivity (11.211 ) under conditions (11.231) and (11.241 ) on a(x): we show that the Wiener condition 
(11.231) is necessary and sufficient for the positivity (11.211) under assumption (11.241) (Theorem [831). 

IV. An asymptotics of the ground state as e —)• 0. 

V. An example of negative energy when the condition (11.241) breaks down while the Wiener condition (11.231) 
holds (Lemma flO.il) . 

VI. Spectral resolution of nonself adjoint Hamilton generators and stability of the linearized dynamics. 

Remark 1.1. The condition ( li.24l) cancels a negative contribution to the energy, which is due to the electrostatic 
instability (’’Earnshaw Theorem” ( 1271/ . see Remark ll0.2\l . 

Let us comment on previous results in these directions. 

The crystal ground state for the Hartree-Lock equations was constructed by Catto, Le Bris, and Lions ||5l|6l. 
Lor the Thomas-Lermie model similar results were obtained in 0|. 

The corresponding ground state in the Schi'ddinger-Poisson model was constructed in llT4ll . The stability for 
the linearized dynamics was not established previously in any model. 

In O, Gances and Stoltz have established the well-posedness for local perturbations of the ground state density 
matrix in an infinite crystal for the reduced Hartree-Lock model of crystal in the random phase approximation 
with the Coulomb potential w(x — y) = l/|x —y|. However, the space-periodic nuclear potential in the equation 
||3l (3)] does not depend on time that corresponds to the fixed ions’s posifions. Thus fhe back reacfion of the 
electrons onto the nuclei is neglected. 

The nonlinear Hartree-Lock dynamics for compact perturbations of the ground state without the random phase 
approximation is not studied yet, see the discussion in ifTTl and in the introductions of ||2l|3l. 

The paper ||2l deals with random reduced HL model of crystal when the ions charge density and the electron 
density matrix are random processes, and the action of the lattice translations on the probability space is er- 
godic. The authors obtain suitable generalizations of the Hoffmann-Ostenhof and Lieb-Thirring inequalities 
for ergodic density matrices, and construt random potential which is a solution to the Poisson equation with the 
corresponding stationary stochastic charge density. The main result is the coincidence of this model with the 
thermodynamic limit in the case of the short range Yukawa interaction. 

In ll2n . Lewin and Sabin established the well-posedness for the reduced von Neumann equation with density 
matrices of infinite trace for pair-wise interaction potentials w G The authors also proved the asymp¬ 

totic stability of the ground state for 2D crystals ll22l . Nevertheless, the case of the Coulomb potential in 3D 
remains open. 

The spectral theory of the Schrddinger operators with space-periodic potentials is well developed, see ll24l and 
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the references therein. The scattering theory for short-range and long-range perturbations of such operators was 
constructed in ifTOlfTTI . 

The plan of our paper is the following. In Section 2 we recall our result ifldll on the existence of a ground 
state, and in Section 3 we establish small charge asymptotics of the ground state. In Sections 4-6 we study 
the Hamilton structure of the linearized dynamics and establish the energy bound from below. In Section 7 we 
calculate the generator of the linearized dynamics in the Fourier-Bloch representation. In Section 8 we prove 
the positivity of energy. In Section 9 we apply this positivity to the stability of the linearized dynamics. Finally, 
in Section 10 we construct examples of negative energy. Appendices concern some technical calculations. 

Acknowledgments The authors are grateful to Herbert Spohn for discussions and remarks. 


2 Space-periodic ground state 

Let us recall the results of ifldll on the existence of the ground state (11.71) . The Poisson equation (11.91) for the 
F-periodic potential d>° implies the neutrality of the periodic cell = M^/F: 

[ p^{x)dx = 0, (2.1) 

which is equivalent to the normalization condition 

f \Y^\x)\^dx = Z ( 2 . 2 ) 

by (11.11) . We assume that Z > 0, since otherwise the theory is trivial. The existence of the ground state (11.71) is 
proved in |[T4]| under the condition 

Op<,r{x) :=Y,<y{x-n) eL^{T^). (2.3) 

n 

The ion position e can be chosen arbitrary, and we will set = 0. 


2.1 Minimization of energy per cell 

The wave function is constructed as a minimal point of the energy per cell 

= \ + p{^)Gp^vp{x)]dx, 

where 

p{x) :=Op,r{x)-e\Y{x)\^, 
while the operator Gper := —-^per is defined by 


Gper^ix) 


y -ilTCmx 


(p{m) = [ (p{x)dx. 


(2.4) 

(2.5) 


( 2 . 6 ) 


More precisely. 


where ^ denotes the manifold 


[/(/) = minUiy/), 


■= {y/e H\T^) : [ \Y{x)\^dx = Z}. 


(2.7) 


( 2 . 8 ) 
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2.2 Smoothness of the ground state 

The results |[T4l imply that there exists a ground state with E Hence E and the 

equation (11.81) implies that 


In ofher words. 

^0 E//4(y.3^ 

)CC^{T^). 


(2.9) 


I {^nf 


(m) := (1 + 

(2.10) 


3 Small-charge asymptotics of the ground state 

We will need below the asymptotics as e —)• 0 of the ground state (11.71 ) corresponding to a one-parametric family 
of ion densities 

a{x)=eiJ.{x) (3.1) 

with some fixed function jj. E We assume that 

jUperW := Y. At(x-n) EL2(r3) (3.2) 

in accordance with (12.31) . Now the energy (12.41) reads 

U{v) = ^^J|Vi/r(x)|He^v(x)GperV(x)]r/x, v(x) := ftperW - \vix)\'^- (3-3) 

Denote by \lf^J,COg the family of ground states with the parameter e E (0,1]. The energy (13.31) is obviously 
bounded uniformly in e E (0,1] for any fixed i/r E Hence, the energy of the minimizers is also bounded 
uniformly in e E (0,1]. In particular, the family i/r® is bounded in (T^), 

||V^“||/r.(r3)<C, ^.E(0,1]. (3.4) 

On the other hand, 

J v”(x)GperV°(x)dx < C, v°(x) := Atpei-(x) - \Yeix)\^- (3-5) 

This estimate is due to the uniform bound 

||v°||i2(r3)<C, eE(0,l] (3.6) 

which holds by (13.21) and (13.41) . Further, the equation (11.91) reads 

-A<I)°(x) =ev°(x). (3.7) 

We will choose the solution d>° = eCpe^Vg, where the operator Gper is defined by (12.61) . The definifion (12.61) 
implies (he bound 

< ^’l|v°L2(y3) < Ce, e E (0,1] (3.8) 

by (13.61) . 

Lemma 3.1. Let condition ( 13.2D hold. Then for sufficiently small e > 0, 

//^=-^A-ed>°(x)-m°>0, (3.9) 

and the ground state admits the following asymptotics as e ^0: 

w° = ^(e2), (3.10) 

¥e{x) = 7e + Xe{x), \Ye\'^ = Z + {e^), \\x ,\\^2 {e^). (3.11) 
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Proof i) Equation (11.81) reads 


(3.12) 


Hence, 


Ve ix) = -2^¥e (x) - (^) {x). 


= <z= -(VvA^,Vt/r«)^3 


which implies the uniform bound 


|w°|<c<oo, ee(0,l] 


(3.13) 

(3.14) 


by (12.21 ). (13.41 ) and (13.81 ). Moreover, (13.121 ) and (13.81 ) suggest that g)° is close to an eigenvalue of — ^A: 

m°«|27r)tp (3.15) 

with some k^l?. Indeed, (13.121) can be rewritten as 


and hence. 


(-|27rm|2 - m°)v>-°(m) = r,(m), r, := 


L (^|27rm|2-mO)2|,)fO(m)|2 = ^(e^), 


(3.16) 


(3.17) 


since ||r(,||^2(7'3) = by (13.81) . Denote by Xe the value of |27rmp corresponding to the minimal magnitude 
of I i|27rnip — (0g\. Then (13.171) implies that 


\27lm\'^^Xe 


(3.18) 


since the set of possible values of ^|27rmp — ft)° is discrete and possible values of ft)° are bounded by (13.141) . 
Moreover, (13.171 ) can be rewritten as 


{^Xe-(0°)^Z+ £ {^\27lm\^-(0°f-{^Xe-(0°f]\¥^{m)\^ = ff{e'^) 




Since 


due to the normalization (12.21) . Hence, 


L \¥eH\^=Z 

^Xe-CO°l = ^(e^), 


since the sum in (13.191) is nonnegative. Let us show that (13.191) also implies that 


^ {\2nm\^-Xef\¥^{m)\^ = ff{e'^). 

\2Km\'^y^Xe 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


First, (13.191 ) gives that 


£ (|2;rm|2-A,)(^|2;rm|2 + lA,-2mO)|i/r; 

\27tm\^^Xe ^ ^ 




However, 2(Og = Xe + G{e^^ by (3.21). Hence, 


^ {\2nm\^-Xe)i\2nm\^-Xe + G{e^))\\f/^{m)f = 

\2Km\'^^k^ 


1 




















Now (13.221) follows from (13.181) since Xe is bounded for small c > 0 by (13.211 ) and (13.141 ). 
ii) Now let us prove that Xe = 0 for small e > 0. Indeed, the energy of the ground state reads 

U{Ve) = i L |27rm|2|i/f0(m)|2 + ^(e2) (3.23) 

^ m€l? 

by (13.31) and (13.51) . On the other hand, (13.221) implies 

^|27rm|^|t//°(m)p = AeZ+ ^ i\2nm\^-Xe)\YKm)\^ = ^eZ+ff{e‘^). (3.24) 

m [2;rmp^Af 

Substituting (13.241 ) into (13.231) . we obtain 

= \XeZ+ff{e^), Xe > 0. (3.25) 

On the other hand, taking t/r(.r) = -\/Z> we ensure that the energy minimum (12.71) does not exceed Hence, 

(13.251) implies that Xe = 0 for small e> 0, since the set of all possible values of XgZ is discrete. Therefore, (13.101) 
holds by (13.211) . 

hi) Now we can prove the asymptotics (13.111) . Namely, the first identity holds if we set 

ye = WeiO), Xe{x) = ^ . (3.26) 

m^O 

Then the second asymptotics of (13.111) holds by (13.201) and (13.181) with Xe = 0. The last asymptotics of (13.111 ) 
holds since 

£ |27rmHvf0(m)|2 = ^(c4) (3.27) 

m^O 

due to (13.221) with Xg = 0. Finally, (13.81) and (13.101) with small c > 0 imply that the lowest eigenvalue of the 
Schrbdinger operator Hg in L?{T^) is close to zero. Hence, its zero eigenvalue is exactly the lowest eigenvalue, 
since the spectrum of this operator is discrete. Therefore, the nonnegativity (13.91 ) is proved for small e > 0. ■ 

4 Linearized dynamics 

Let us consider the linearized system (11.131) . We recall that G := — The meaning of the terms with G will 
be adjusted below, see Lemma [531 We assume further that (12.3!) holds, and additionally, 

(x)2aGL2(M^), (A-l)aGLi(M3). (4.1) 

For f{x) G C^(M^) the Fourier transform is defined by 

m = 7^l m)= f e‘^^f{x)dx, ^GM^ (4.2) 

{2ny 7r3 

The conditions (14.11) imply that 

(A — l)a G L^(M^), (i^)^a(i^) < const. (4.3) 

We consider the case when the ground state Y^{x) can be taken to be a real function. Then (IL13l) - (IL15i) imply 
that the operator-matrix A is given by (11.171 ) where S denotes the operator with the “matrix” 


S{x,n) := eY^{x)GV(j{x — n) : n G 1?, x G 


( 4 . 4 ) 













































Finally, T is the real matrix with entries 

T{n,n) := — (GV ®'S/o{x — n'),o{x — n)) + ( 4 >°,V ®Vo)dnn' = T\{n — n) + T2{n — n). (4.5) 

The operators Gyr® : —)• L^(M^) and S : ■ = Z|(Z^) (8)C^ —)• L^{M?) are not bounded due to the “infrared 

divergence”, see Remark [54l In the next section, we will construct a dense domain for all these operators. 

On the other hand, the corresponding operators T\ and T2 are bounded by the following lemma. Denote by 
n the primitive cell 

Yl\= {{x\,X2,X2,) '-0 <Xk k = 1,2,3}. (4.6) 

Let us define the Fourier transform on Z 3 as 

(2(e) = £ e'”®(2(u), a.e. e G n*; Q{n) = / e-‘"^Qid)dG, n G Z^, (4.7) 

|n*l 


where IT* = 27111 denotes the primitive cell of the lattice F* and the series converges in 

Lemma 4.1. The operators Ti and T 2 are bounded in Z 3 under condition ( I4.il) . 

Proof The first operator Ti reads as the convolution: TiQ{n) = l}ri(n — n')Q{n'), where 

Ti{n) =—(V iS>G'Vo{x),a{x — n)). (4.8) 

In the Fourier transform (14.71) . the convolution operator T\ becomes the multiplication, 

T\Qie) = fiie)Q{e), a.e.een*\r*. (4.9) 

By the Parseval identity, it suffices fo check fhaf fhe “symbol” Ti(P) is a bounded function. This follows by 
direcf calculafion from (14.51) . Firsf, we apply fhe Parseval idenfify: 


Me) = -£.-e(V®GVa(v),a(v-n)) = ^£.'«^^a(^),a(^).-^) 






L 




- ^=2nm—d 


, eGn*\r* (4.io) 


(27r)3' 

since fhe sum over n equals |n*| ^5(0 +— Inm) by fhe Poisson summation formula ifT^ . Finally, |^(i^)| < 

m 

by (14.31) . Hence, 


|Zi(0)| <Ci£|a(27rm-e)a(27rm-0)| < C 2 £(m)^'* < oo. 

m m 


(4.11) 


ii) Finally, 

7^(0) = f 2 ( 2 ( 0 ), 0Gn*, (4.12) 

where 

f 2 = (d>‘’(x),V®Va(v)) (4.13) 

by (11.9b . The expression is finite by (14.1b . since <I>‘^ G is a bounded periodic function. ■ 
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5 The Hamilton structure and the domain 


To construct solutions of the system (11.171) . we need to diagonalize its generator A. The main problem is 
that this generator is neither selfadjoint and even not symmetric, so we cannot apply the von Neumann spectral 
theorem. We will solve this problem by applying our spectral theory of Hamilton operators with positive energy 
ifTSlfT^ to the Bloch representation of A. 

In this section we study the domain of the generator A. Denote 

r (5.1) 


It is easy to check that the Hamilton representation (11.181) formally holds with the symplectic matrix 


J = 


I 0 

2 

0 
0 


1 

2 

0 

0 

0 


0 

0 

0 

-1 


o\ 
0 
1 

0 / 


(5.2) 


Definition 5.1. i) := UooA^e, where is the space of functions (p G .^^(M^), whose Fourier transforms 
vanish in the e-neighborhood of the lattice T*, 

ii) Ic = UKgNZc(-^). where lc{R) '■= {G G Z 3 : Q{n) = 0 for \n\ > R}. 

iii) ^:={Y = {^'i,^'2,Q,P)(^^:^u^2^y+, Q&lc, P^lc}. 


Obviously, ^ is dense in 

Theorem 5.2. Let conditions ( I4.il) hold. Then B is a symmetric operator on the domain C X. 

Proof Formally the matrix (11.181) is symmetric. The following lemma implies that B is defined on Lf. 
Lemma 5.3. i) G L^(ffi^) and S*(p G l^for tp G 

ii) SQ£L^{m?)forQ£F. 


Proof i) First, note that 




Further, = (27r)^LmGZ3 —Inm). Respectively, 


(5.3) 


= (271)3 £ t/f‘’(m)(p((§-27rm) =0, \^\< £ (5.4) 

if (p G A^g with some e > 0. Moreover, £ L^(M^), since t/r°(p G L^(M^). Hence, (p belongs to the domain 

of G\f/^ and of 

Now consider S*(p. Applying (14.41) . the Parseval identity and (15.41) we get for cp G 


[S*cp]{n) 


e 


\l/^{x)(p{x)GVG{x — n)dx = e{\l/^{x)(p{x),GVa{x — n)) 


le 


(271); 


i 


^\>e Isl 

Here * (p]{^){^)^ G L^(M^) for all a by (12.101) . since 9 G A^(M^). Moreover, d^d 

by (14.31) . Hence, integrating by parts twice, and taking into account (15.41) . we obtain 


|[S>](n)|<C(n) 


(5.5) 

G L^(M^) for |a| < 2 

(5.6) 
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which implies that S*(p G 

ii) Let us check that SQ G L^(M^) for Q G Ic- The Fourier transform of SQ reads as 

n n 

= e{2n)^ [Y,¥^(>n)d{ri -2nm)GVa(^ -ri)Y^e‘'^^^^^^Q(n)dri 

m n 

= e(2;r)3^t/fO(m)G^(i^-2;rm)(2(i§-2;rm). (5.7) 

m 

Hence, the Parseval identity gives that 

||5e|L2(R3) =C||^|L2 (r 3) <Ci£|vrO(m)|||GV^(^)(2(^)||^3(K3) (5.8) 

m 

It remains to note that the sum over m is finite by (12.101) because 

||GV^G|li2(K3) = f j^\a{^)m\^d^<C{Q) p-^l^d^ (5.9) 

since the function Q{^) is bounded for Q G Ic- Finally, the last integral is finite by (14.31) . ■ 

This lemma implies that BY G ^ for Y G The symmetry of B on is evident from (11.181) . Theorem l5.2l is 
proved. ■ 

Remark 5.4. Tlw infrared singularity at ^ = 0 of the integrands ( 15.51) . ( 15.5D and ( 15.9D demonstrates that all 
operators Gy/^ : —)■ L^(M^), 5 : L^(M^) and S* : L^(M^) —)• l^ are unbounded. 

Corollary 5.5. The proof of Theorem 15.21 shows that the operator A is defined on 3!, as well as the ’’formal 
adjoint” A*, which is defined by the identity 

{AY,,Y2) = {YuA*Y2), YuY2G3. (5.10) 


6 Factorization of energy and bound from below 

The equation (11.171 ) is formally a Hamiltonian system with Hamilton functional j{Y,BY). Next theorem means 
the stability property of the linearized crystal. 

Theorem 6.1. Let conditions ( 14.71) hold. Then the operator B on the domain 3) is bounded from below: 

{Y,BY)>-C\\Yf^, Yg3. (6.1) 

Proof For Y = (fYi,'Y2,Q,P) ^ ^ the quadratic form reads (with the notations (I4.4I) - (I4.5I) ) 

(T,BT) = 2Y^{^>j,H^^>j)AAe^{^p^^>uG^p^^>,)+2[{^>l,SQ) + {Q,S*'Yl)] + {QJ^Q) 

j 

+ {Q,T2Q) + {P,M-^P)- ( 6 . 2 ) 

Here the first sum is bounded from below, the operator T 2 is bounded in Z| by Lemma |4~T1 while the operator 
is positive. Our basic observation is that 

i3('Pi,G) := 4 c2(/'Pi,Gva‘’'Pi) +2[('Pi,50 + ((2,5*'Pi)] + {Q,TQ) > 0. (6.3) 
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Indeed, the operators factorize as follows: 


eV‘’G/=/7, s=rg, n=g*g, (6.4) 

where 

f:=eVG\lf^\ g{x,n) = 'VVGo{x — n). (6.5) 

Then the quadratic form (16.31) becomes the ’’perfect square” 

l5{^',Q) = {2f^i+gQ,2f^i+gQ)>0. m (6.6) 

Corollary 6.2. The operator B admits selfadjoint extensions by the Friedrichs Theorem / I23I/ . 


7 Generator in the Fourier-Bloch transform 

We reduce the operators A = JB and K by the Fourier-Bloch-Gelfand-Zak transform |[7ll26l. 


7.1 The discrete Fourier transform 

Let us consider a vector Y = (*Ti,* 1 ^ 2 ,2,P) G and denote 

Y{n) = {'¥i{n,-),'¥ 2 {n,-),Q{n),P{n)) , n^I?, (7.1) 


where 


'Ty(n,y) = 


_ / 'T,(?i+y), a.e. y G H, 


0 , y^n. 

Obviously, Y [n) with different n G are orthogonal vectors in ^, and 

Y = l,Y{n), 


(7.2) 


(7.3) 


where the sum converges in The norms in ^ and y can be represented as 


I rWII 

nCil? 


2 

.r(n)i 


\Y\?r= I \\Y{n)\\ 

nCil? 


2 

r(n)> 


(7.4) 


where 

,r(n) :=L2(^)©L2(^)©c3©c^ r(n) :=//'(n)©//'(n)©c3©c^ (7.5) 

Further, the ground state (11.71) is invariant with respect to translations of the lattice F, and hence the operator A 
commutes with these translations. Namely, (14.41) implies that 

S{x,n) = S{x — n,Q), (7.6) 

since V^®(v) is a F-periodic function. Similarly, (14.51) implies that T commutes with translations of F. Hence, A 
can be reduced by the discrete Fourier transform. Namely, applying the Fourier transform to the function 
T(-) from (TtH) . we obtain 


f(0)=F„^eZ(n):= £ c'”7(n) = ('Fi(e,-),'F2(0,-),e(e),P(0)), a.e. 6 G (7.7) 

where 

+ a.e. 0GM^ a.e. y G (7.8) 
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The function Y(d) is F*-periodic in 6. The series (17.71) converges in L^(n*,^(n)), since the series (17.31) 
converges in The inversion formula is given by 


Y{n) = \n*\-'^ [ e-’"^Y{d)dd 
Jn* 

(cf. (14.71) ). The Parseval-Plancherel identity gives 


lly ||2 _ ITT* l|y 11 ^ 

IP \\r — I IP \\v-{u*.r(n))i 

The functions ^j{0,y) are F-quasiperiodic iny; i.e., 


I^l|jr — |n I \\^\\L^(n*,s:{u)y 


T*| —1 llvl |2 




(7.9) 


(7.10) 


(7.11) 


7.2 Generator in the discrete Fourier transform 

Let us consider Y £ ^ and calculate the Fourier transform (17.71) for AY. Using (14.51) . (15.51) . (17.61) . and taking 
into account the F-periodicity of d>‘’(.r) and we obtain that 


^{d)=A{d)Y{d), a.e. 0eR^\F*, 
where A(0) is a F*-periodic operator function, 

/ 0 0 0 \ 

-H^-2e^Y^G{d)Y^ 0 S{d) 0 


(7.12) 


A(0) = 


V 


0 

- 2 s*{e) 


0 0 M-^ 

0 -f{G) 0 y 


(7.13) 


by (11.171) and (11.181) . Here 


j(2nm-e)y 




a.e. G£R^\r* 


This expression is well-defined for q){x) = i/r°(x)'Pi(x) with 'Ll £ since 

(p{G,m) = ^(2nm —G) =0 for |27rm —0|<£ 

according to (15.41) . 

Lemma 7.1, The operator S{G) acts as follows: 

S{G)Q{G) = S{G)Q{G), where 5(0) = ei//°G(0)Vd(e,y). 
Proof. For x = y + n equations (11.111) and (14.41) imply 

SQ{y + n) = e'^^{y + n)'Y^GVo^{m,y + n)Q{m) 

m 

= e\ff^{y)Y,GVa{y + n — m)Q{m) 

m 

due to the F-periodicity of i/r®. Applying the Fourier transform (17.71) . we obtain (17.161) . 


(7.14) 


(7.15) 


(7.16) 


□ 


Furthermore, S*{G) in (17.131) is the corresponding adjoint operator, and t{G) is the operator matrix ex¬ 
pressed by (14.101) . Note that S(0), S*{G) and T{G) are finite dimensional operators. 
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7.3 Generator in the Bloch transform 


Definition 7.2. The Bloch transform ofY is defined as 

f(e) = [^F](e) := ^{0)Y{e) -.= {^yiO,y),%id,y),Q{d),P{0)), a.e. d € M', 

where j) = M(0)'Fy := e^^^’^j{6,y) are T-periodic functions in y S 
Now the Parseval-Plancherel identities (17.101) read 


|y||2 _ |TT*l^*llyll^ 

P Hr — I IP llL2(n*,r(r3))) 


Y\\1r = \^ 


*1—1 I 


llL2(n*,.r(7'3))- 


Hence, ^ \ ^ L^{Yl*, ^{T^)) is, the isomorphism. The inversion is given by 

Y{n) = \n*\-^ [ e-‘"^^{-e)Y{d)de, n£l?. 

in* 

Finally, the above calculations can be summarised as follows: (17.121) implies that for Y G 

^( 0 ) =A(e)T( 0 ), a.e. eGn*\r*. 

Here 

A(e)=^(e)A(e)^(-e) = 


where 



/ 

0 

H^{d) 

0 

0 \ 

6 ) = 



0 

SiG) 

0 


0 

0 

0 



1 

-2S*{d) 

0 

-f( 6 ) 

0 / 


Remark 7.3. The operators (5(6) : lf{T^) — )■ H^iT^) are bounded for 6 G H* \r*. 

Lemma 7.4. Let the condition ( 17.271) hold. Then the operator A{d) admits the representation 

Aid) = jB{d), eGn*\r*, 
where B{d) is the self adjoint operator ( 17.791) in ^ {T^) with the domain 

b := 77^(r^) ©77^(r^) © © cl 


(7.17) 


(7.18) 


(7.19) 


(7.20) 


(7.21) 


5(0) 

:= M(0)5(0) =et/r‘><5(0)Va‘’(0), 

(7.22) 

^°( 0 ) 

:= M{e)H^M{-e) = -^(V + /0)2 - e^^ix) - ml 

(7.23) 

(5(0) 

:= M( 0 )G( 0 )M(-e) = (/V- 0 )-l 

(7.24) 


(7.25) 


(7.26) 


Proof The representation (17.251) follows from (11.181) and (11.181) . The operator B{6) is symmetric on the domain 
b. Moreover, operators in (11.191 ) are all bounded, except for which is selfadjoint in L?{T^) with the 

domain H^{T^). Hence, B{d) is also selfadjoint on the domain D. ■ 
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8 The positivity of energy 


Here we prove the positivity (11.211) for the linearized dynamics (11.171) under conditions (11.231) and (11.241) . It is 
easy to construct the corresponding examples of densities a(x). 

Example 8.1. (17.231) holds for o{x) € L^(M^) if 


d{^) / 0 , a.e. ^ e 


( 8 . 1 ) 


0 ■ ^ 

2 sm — 

~ O e2 

Example 8.2. Let us define the function f{x) by its Fourier transform f{^) := —^ 

a{x) :=eZf{xi)f{x 2 )f{x 3 ), xgR^. (8.2) 

Then o{x) is the smooth function satisfying the Wiener condition ( 17.231 ). as well as ( 17.241) and ( 17.71) . and 

|a(x)| <C(a)e-"W, (8.3) 

for any a> 0 by the Paley—Wiener theorem. 

The matrix (11.231 ) is a continuous function of 0 G IT* \ F*. Let us denote 

n; := {0 Gn*\r* :i( 0 ) > 0 }. (8.4) 

Then the Wiener condition (11.231 ) means that |n'*j_| = |n*|. In the rest of this paper we assume condition (11.241 ) 
and consider the linearized dynamics (11.171) corresponding to a real minimizer of energy per cell. In Appendix 
B we show that the real minimizer exists and is unique. 

Theorem 8.3. Let conditions ( 14. 7 D . and ( 17.241 ) hold. Then the Wiener condition ( 17.231 ) is necessary and suffi¬ 
cient for the positivity ( 17.271) of the generator corresponding to the real minimizer of energy per cell. 

Proof i) First, let us check that the Wiener condition (11.231 ) is necessary. Namely, let us consider the inequality 
(11.211) for To = (0)0)G)^) G 'V{T^)'. (11.191 ) and (11.211) imply that 

S{e,YQ) = Qf{e)Q + PM-^P>K{e)[\Q\^ + \PW a.e. 0Gn*\r*. (8.5) 

Formula (14.131) implies that 72 = 0 by (IB.51) . Hence, 

f(0) = fi(0) = r(0), 0Gn*\r* (8.6) 

by (14.101 ). Therefore, (18.51 ) becomes 

^(0,To) = Gr(0)G + PM-ip > x(0)[|(2|2 + |P|2]. (8.7) 

Hence, the condition (11.231) is necessary for the positivity (11.211) . 

ii) It remains to show that the Wiener condition (11.231) together with (11.241) is sufficient for the positivity (11.211) . 
Let us translate the calculations (I6.2I) - (I6.5I) into the Fourier-Bloch transform. The operators (16.51) commute 
with the F-translations, and therefore 

^>V‘’G( 0 )t/rO=f( 0 )/( 0 ), (8 8 ) 

where /(0) := e-\/ G(0)t/r‘’ and g(0) = Y^(5(0)Va(-, 0). Hence, (11.191 ) implies that 
S{e,Y)-.= {Y,B{e)Y)T,=b{e,W^,Q)+2{W2,H\e)W2)T^+PM-^P, F = ('Fi,'F2)G)P)Gr(p3), (8.9) 
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where 


b{e,^>,,Q) ■=l{^>uH\^)^>,)T^ + {i~m^, + ~g{d)Q, 2/(e)'Fi+i(e)0r3. (s.io) 

Let us note that H^{d) = —-(V + /6)^ by (IB.5I) . Hence, the eigenvalues of ^*^(0) equal to -|27rm — 6\^ where 
m^7?. Therefore, H^{d) is positive definite: 

('Ti,^0(e)'Ti)>^d2(e)||'Ti||^,(^3), 0 Gn*\r, (s.ii) 

where d{d) := dist (0,r*). Hence, it remains to prove the following proposition. 

Proposition 8.4. Under conditions of Theorem \8.3\ 

M0,'Pi,G)>£(0)[||'Pi||^i(r3) + |Gp], Sen;, (8.12) 

where e{6) > 0 . 

Proof Let us denote a := (*Pi ,^°(0)*Pi)7'3, and 

^n:=(2/(0)'Pi,2/(e)'Pi)r3, ^12 := (2/(e)'Pi,g(0)0r3, fe := (|(0)0g(0)0r3. (8.13) 

Then we can write the quadratic form (18.101) as 

b = 2a + l5, j8 :=i3ii+2Rej8i2 + i 822 . (8.14) 


The positivity (18.111) implies that 
where 5{d) > 0. Hence, 


a>5(0)j3ii, 0Gn*\r, 


(8.15) 


Z7>a + (l + 5(0))jSii+2Rej8i2 + i 822 , 0 Gn*\r. (8.16) 


On the other hand, the Cauchy-Schwarz inequality implies that 

1^121 <^;fi82f <^[7^11+^)822] 

z 7 

for any 7 > 0. Hence, 

b>a + {l + 6(9) - 7)jSii + (1 - ^)p22, 9 G H* \r. 

Therefore, choosing !</< 1 + 5(6), we obtain (18.121) from (18.1 II) since 

(822 = Gri(0)(2 = 1(6)102 


by (18.81) and (18.61) . 


(8.17) 


(8.18) 


(8.19) 


9 Weak solutions and linear stability 

Weak solutions are introduced and the linear stability is proved. 
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9.1 Weak solutions 


We will consider solutions to (11.171) in the sense of distributions. Let us recall that A*V G ^ for L e ^ by 
Corollary 15.51 

Definition 9.1. Y{t) S C(M, is a weak solution to if 

- J{Y{t),(p{t)V)dt = J{Y{t),(p{t)A*V)dt, (peCo(M), (9.1) 

Let us translate this definition into the Fourier-Bloch transform: by the Parseval-Plancherel identity 


LJn* 


{Y{G,t),(pit)V{d))r^dd 


dt = 


/n» 


{Y{ef)Mt)A*{e)vmT^de 


dt 


(9.2) 


Respectively, (19.11 ) is equivalent to the identity 


,t),(pit)V)T3dt = 


,t),^{t)A*(G)V)T3dt, (OGCq 


Ve^iT^), a.e. 0Gn*\r, (9.3) 


where ^(r^) := C“(r3) ©C“(r3) © © Cl In other words, 

f{G,t)=A{G)Y{Gf), a.e. eGn*\r (9.4) 


in the sense of vector-valued distributions. 


9.2 Linear stability 

The equation (19.41 ) is equivalent to 

f{Gf)=JB{G)Y{Gf), fGM, a.e. 0Gn*\r. (9.5) 

We reduce it, using (11.211) . to an equation with a selfadjoint generator by our methods ifTSl [T^ which is an 
infinite-dimensional version of some Gohberg and Krein ideas from the theory of parametric resonance ^V2\ 
Chap. VI]. We reproduce some details of |[T5l for the convenience of the reader. Namely, let us denote 

A(0) =bV 2(0) >0, PGn;. (9.6) 

This is a selfadjoint operator with the domain y{T^), that follows by the interpolation arguments, and the range 
Ran A(0) = {T^). Its inverse is bounded in (T^) by (11.211) . and 

||A-i(0)Z||^(r3)<^i=||Z||^(r3), ZGl%r{T^), GGUf. (9.7) 

V 

Let us set Z{G,t) := A{G)Y{G,t), and now equation (19.51) implies that 

Z{G,t) = -iK{G)Z{G,t), teR, a.e. 0 G n+ (9.8) 

in the sense of vector-valued distributions, where K{G) = iA{G)JA{G). 

Lemma 9.2. (Lemma 2.1 of ifTSl l K{G) is a selfadjoint operator in ^{T^) with a dense domain D{K{G)) C 
f'{T^) for every G G HI. 
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Proof The operator K{d) is injective. On the other hand, Ran A(6) = ^(T^), and J : X(T^) —)• X(T^) is a 
bounded invertible operator. Hence, Ran.^(0) = X{T^). Consider the inverse operator 

R(e) ■=K-\e)=iA-\e)JA-\d). (9.9) 

It is selfadjoint since D{R{d)) = Ran A'(6) = X(T^) and R{d) is bounded and symmetric. Finally, R{d) is 
injective, and hence, .^(0)=R^*(0)isa densely defined selfadjoint operator by Theorem 13.11 (b) of |[25l : 

K*{d)=K{d) , D{k{e)) =RmR{d) C RanA-i(e) C y{T^) 

bv (lO). ■ 

This lemma implies that each weak solution to (19.81 ) is given by 

z(e,t) GCfo(M,s'(r3)), a.e. een; ( 9 . 10 ) 

for Z(0,O) G X (T^). Hence, we obtain the well posedness of the Cauchy problem for equation (19.51 ). 

Theorem 9.3. Let all conditions ofTheorem \8.3\ hold and 6 G H’^. Then for every initial state Y{d,0) G iC 
there exists a unique weak solution Y{d,t) G Cfo(M, k{T^)) to equation ( 19.51) . and 

(A(e)T(e,t),A(0)T(0,f))7.3 = const, t G M. (9.11) 

Proof Z(0,O) := A{6)Y (0,0) G X(T^) since Y (0,0) G y{T^). Hence, (19.101 ) and (19.71) imply that 

Y{ef)=A-\d)e-'^^^^*Z{d,0)eCb{'R,y{T^)). (9.12) 

Finally, (19.1 II ) holds since is the unitary group in X(T^), and hence 

(A(0)T(0,t),A(0)T(0,t))r3 = (Z(0,t),Z(0,t))r3 = const, f G M. ■ 

Now we apply this theory to equation (11.171) . Let us note that A(0)T(0) G L^(nf,X (T^)) for Y G see 
Definition l5.ll 

Definition 9.4. The Hilbert space "W is the completion of in the norm 

||F||^:=||A(0)T(0)||i2(n;,^(r3)) (9.13) 

Formally, ||T||^ = {Y,BY)^^^. The Fourier-Bloch transform (17.171 ) extends to the isomorphism 

^■.W^W-.= {Y{-)glUuI,X{T^)) : ||A(0)f(0)|b(n*,^(r3)) <-}. (9.14) 

Finally, let us extend definition of weak solutions to Y(t) G Cb(R,^) by the identity (19.51) in the sense of 

vector-valued distributions (19.31 ). Then Theorem |93] implies the following corollary. 

Corollary 9.5. Let ail conditions of Theorem 18.51 hold. Then for every initial state Y{Q)gW there exists a 
unique weak solution Y(t) G Cfo(M, W) to equation ( 17.771) . and the energy norm is conserved: 

||T(f) 11^ = const, t G M. (9.15) 

The solution is given by the formula ( 19. 7 2D .' 

T(f) = ^-iA-'(0)e-'^(®)'Z(0,O) GCi(M,^(r3)). (9.16) 

This means that the linearized dynamics (I1.17D is stable: global solutions exist for all initial states of finite 
energy, and the norm is constant in time. 
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10 Examples of negative energy 


We show that the positivity (11.211) can fail if the condition (11.241) breaks down even when the Wiener condition 
(11.231) holds. Namely, for Yq = (0,0,2,0) G (T^) we have 


c^iejo) = Qf (6) Q 


( 10 . 1 ) 


by (IQ) . 


Lemma 10.1. There exist functions }x{x) such that the positivity ( li.2il) fails for cj(.r) from ( I3.il) with small 
e > 0 while ( 14.il ) and the Wiener condition ( li.23D hold. 


Proof It suffices to construct an example of o{x) which provides 

Gf(eo)G<0 

for some 6o G n* \ F* and The representation (14.101) can be written as 


Ue) = e^Y. 




-\m)Y 


. ^= 271111-9 


, een*\r* 


Similarly, (14.131) can be written in the Fourier representation as 


T2 = -e" 




(2;r)3' 


( 10 . 2 ) 


(10.3) 


(10.4) 


with v° := itper(-^) — \Ve(.^)\^ according to (13.51) . The asymptotics (13.111) of the ground state \f/^J{x) implies 
v°(§) = Atper(^)-|7.P(27r)'5(^)-.~(§) = Atper(^)-Z(27r)35(^)-.~(§), 
since = Z by ( 13.111 ). Here s{x) = YeXei^) + 7 eXe{x) + |ZeWP> and 

lk||L2(7.3) < 

by (13.111 ). Further, (13.21 ) gives 

Aper(^) = lAi(^)^'"^ = At(^)(27r)3£5(§ -27tm) 

n m 

by the Poisson summation formula ifT^ . Substituting (110.71) into (110.51) we get 

v0(^) = A(^)(27r)3£5(^-27rm)-.~(^) 

rtiy^O 

by (11.11) and (13.11) . Substituting this expression into (14.131) we obtain 

f2 = -e^m) 15 (^ -27im)^rm) 


(10.5) 


( 10 . 6 ) 


(10.7) 


( 10 . 8 ) 


(10.9) 


m^O 


At last, s{x) is a F-periodic function and 


j^,s{x)dx = l^^v^ix)d. 


by (13.71) . Hence, 


'x = 0 


^|5(m)|^ = (^(e^), e —)■ 0 

m^o m 


( 10 . 10 ) 
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by (I10.6I) . Therefore, 


f2 = -e^E 

m^O 






\m)? 


. ^=27Cm 




0 . 


( 10 . 11 ) 


Hence, there exists a 2 € such that 

Qf2Q<0 ( 10 . 12 ) 

for small e > 0 if the condition (IB.21) breaks down. For example, we can take Q = 2nm with m S \ 0 if 
jl(2Km) / 0. Finally, for any Oq 0 T* we can reduce |/t(0)| in all points 0 S Oq + r* keeping it in the points of 
F* to have 

Qfido)Q = Qt,{eo)Q + Qt2Q<0. (10.13) 

At the same time, we can keep (14.11) and the Wiener condition (11.231) to hold. ■ 


Remark 10.2. The operator T 2 corresponds to the last term in the last line of ( 17.731) . This term describes the 
’’virtual repulsion ” of the ion located at n + from the same ion deflected to the point n + q^ + Q{n,t)- This 
means that the negative energy contribution is provided by the electrostatic instability (’’Eamshaw Theorem ” 
1271). 

A Formal linearization at the ground state 

Let us substitute 

V/(x,t) = [v/°(x)+'F(x,t)]e^'®°', q{nd) =q^ + Q{n,t) 
into the nonlinear equations (11.21) . (11.41) with ^{xf) = Gp(x,t). First, (11.31) implies that 

p{x,t) = '^a{x — n — q^ — Q{n,t)) —e\Y^{x) +'F(x,t)p 


and the Taylor expansion/orma/Zy gives 


p{x,t) = ^ o{x — n — q^) — 'Vo{x — n — q^)Q{n,t) + -'V'Vo{x — n — q^)Q{n,t) ^Q{n,t) + ... 

|V^‘’(x)|^ + 2Re(t/r'’(x)'F(x,t)) + |'F(x,t)p = p°(x) +pi(x,i) +p 2 (x,t) + ... 


— e 


(A.l) 


Here p°(x) := a‘’(x) —e|v/°(x)|^ and pj^ are polynomials in 'F(x,t) and Q{t) of degree k. In particular, p\{xf) 
is given by (11.141) . As a result, we obtain the system (11.131) in the linear approximation. 


B Ground states with minimal energy per cell 

Let us consider any ion density o{x) G L^(M^) satisfying (11.241) : 

a(27rm)=0, mGZ^\0. (B.2) 

Let us note that 

a(0) = J a{x)dx = eZ>0 (B.3) 

by (11.11) . Then aper(x) := a(x — n) = eZ since 

ape..(m)= f e‘^^'"^aper{x)dx= f e^^’”’‘a{x)dx = 0, mGZ^\0 (B.4) 

JG Jr^ 


20 























by (IB.2I) . Therefore, the functions 


/(a:) = Vz, 4)‘’(a:) = 0, w‘’ = 0 (B.5) 

give a solution to (I1.8I) - (I1.10I) with zero energy per cell (12.41) . On the other hand, the energy (12.41) is nonnegative. 
Hence, the set of all minimizers of energy per cell consists of \/Z, with <j) £ [0,2;r]. 
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